Many-Impurity Effects in Fourier Transform Scanning Tunneling Spectroscopy 
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Fourier transform scanning tunneling spectroscopy (FTSTS) is a useful technique for extracting 
details of the momentum-resolved electronic band structure from inhomogeneities in the local density 
of states due to disorder-related quasiparticle scattering. To a large extent, current understanding 
of FTSTS is based on models of Friedel oscillations near isolated impurities. Here, a framework for 
understanding many-impurity effects is developed based on a systematic treatment of the variance 
Ap 2 (q, uj) of the Fourier transformed local density of states p(q, w). One important consequence of 
this work is a demonstration that the poor signal-to-noise ratio inherent in p(q, ui) due to random- 
ness in impurity positions can be eliminated by configuration averaging Ap 2 (q, uj). Furthermore, 
we develop a diagrammatic perturbation theory for Ap 2 (q,u) and show that an important bulk 
quantity, the mean- free-path, can be extracted from FTSTS experiments. 
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I. INTRODUCTION 

As a result of technical advances in scanning tunneling 
microscopy (STM) methods, it has become possible to 
measure the quasiparticle spectrum at surfaces of a va- 
riety of materials over large areas with both atomic spa- 
tial resolution and mV energy resolution. While STM is 
a powerful and direct probe of electronic wavcfunctions, 
the interpretation of STM spectra generally requires com- 
parison with a model spectrum calculated using approxi- 
mate many-body techniques (for example, the local den- 
sity approximation). This can be problematic in situa- 
tions where details of the surface are poorly known, or 
where traditional theoretical methods fail, as in the high 
temperature superconductors. The question of how one 
can usefully analyze STM spectral maps without a de- 
tailed quantitative model of the electronic structure is 
therefore paramount. 

Recently, it was shown that details of the quasiparti- 
cle dispersion can be extracted from the spatial Fourier 
transform of STM spectra on disordered surfacesi The 
technique is often referred to as Fourier transform scan- 
ning tunneling spectroscopy (FTSTS) and requires that 
maps of the local density of states (LDOS) p(r, uj) 
at energy u> be made over reasonably large surface 
areas, from which the Fourier transform, p(q, to) = 
J2 r i°( r > w ) exp(— iq • r), can be obtained. In the simplest 
case, that of noninteracting quasiparticles, the LDOS at 
position r and energy uj is 



P(r,w) 



E 



V n (r)\ 2 6(cu-E n 



(1) 



where ^(r) is the n th quasiparticle eigenstate with en- 
ergy E n . For a clean surface, $ n (r) are Bloch states of 
the form 3 - i ] t (r)e flt " r with wavevector k and band index 
j, and p(r,uj) = ^ ■ k \<f>j jk (r)\ 2 S(uj - E jk ) has only the 
spatial periodic modulations of the unit cell. However, in 
the presence of defects, the eigenstates contain mixtures 
of Bloch states with different values of k, and the proba- 
bility density |4' n (r)| 2 exhibits standing wave patterns 
(Friedel oscillations) arising from the interference be- 
tween these different k- values. Near an isolated impurity, 



for example, one can decompose the eigenstates at energy 
uj into initial and final (scattered) states with wavevec- 
tors ki and kf respectively. Then, the Fourier trans- 
formed density of states (FTDOS) is peaked at wavevec- 
tors q = k/ — ki for which k/ and k^ lie on constant 
energy contours of energy uj. This is illustrated in Fig.^ 
in panel (a) the FTDOS is calculated (details given later) 
at the Fermi energy uj = for a single impurity in a 
2D tight-binding metal whose Fermi surface is shown in 
panel (b). The arrows qi and q2 in (b) label transitions 
between states on the Fermi surface which contribute to 
the maxima at qi and q2 in (a). There is, therefore, a 
direct connection between p(q, uj) and the quasiparticle 
energy dispersion ek- The possibility of using FTSTS to 
map out ek was first demonstrated by Sprunger et alX 
and FTSTS has since been used by several groups to 
study metal surfaces)^ semiconductors semimetals,? 
and strongly-correlated superconductorsA&SiiSi Quite re- 
cently, an extension of this technique to study inelas- 
tic scattering in strongly-correlated superconductors has 
been proposed^ 

The usefulness of the FTSTS method, however, is lim- 
ited by the large signal-to-noise ratio which is invariably 
present. This is an intrinsic problem, as we now illus- 
trate, which arises from the randomness of the impu- 
rity positions. We assume, for simplicity, that a ran- 
dom distribution of Ni identical impurities sit at sites 
Rj, i = 1, . . . , Ni and that the impurities do not interfere 
substantially, so that each impurity modifies the LDOS 
by an amount Sp(r — Rj,o;). This term is meant to de- 
scribe the Friedel oscillations in the vicinity of an isolated 
impurity, and its Fourier transform Sp(q, uj) is the quan- 
tity one wishes to study experimentally. For simplicity, 
we will ignore sub- unit-cell modulations of the LDOS due 
to the crystalline potential so that the LDOS is uniform 
in the absence of impurities. Then, with the impurities 



present, p(r,u) = p (uj) + J2 
follows that 



1 Sp(r — Rj,w), and it 



p(q,uj) = [Np (uj) + Ni5p(q,uj)}5^ + 5p(q,uj)^N~f(q), 

(2) 

where 5p(q,uj) — J2r e l<i ' T Sp(r,uj) is the Fourier trans- 
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form of the scattering pattern due to a single impurity, 
N is the number of spatial points included in the Fourier 
transform, and 



/(q) 



s/Nl 



E' 

i=l 



-zq-R,i 



is a random function of q normalized such that |/(q)| ~ 1 
and (5 q ^o = 1 — <5q.o- The form of Eq. J2J was first demon- 
strated in the limit of weak impurities by Capriotti et 
alw^k and then for a dilute distribution of pointlike impu- 
rities of arbitrary strength by Zhu et alm&. We make two 
comments regarding Eq. J2J. First, for a particular dis- 
order configuration, and for q ^ 0, p(q, u) is the product 
of a useful signal <5p(q, lo) and a noise term /(q). In prac- 
tice, this means that the noise is as large as the signal, 
and that it can be nearly impossible to resolve the fine de- 
tails of the FTDOS spectrum which have been predicted 
in numerous calculations. This is particularly true in 
the cuprate superconductors, which are inherently dirty 
materials. This point is illustrated in Fig. where the 
FTDOS is shown for a 2D metal for both a single point 
scatterer in panel (a), and a random configuration of 3% 
impurities in panel (c) (the details of the calculation will 
be given below) . The figure clearly demonstrates the dif- 
ficulty one faces when dealing with disordered systems; 
while the FTDOS pattern is visible in (c), one would have 
difficulty making a quantitative comparison of the struc- 
ture with (a). The second comment is that in nearly all 
calculations to date, £>p(q, lo) has been calcuated in the 
clean limit, ie. for a single impurity. This approximation 
captures many qualitative features of the spectrum, but 
also risks failing when the scattering rate is large, as ap- 
pears to be the case in the cuprate superconductors. It 
is natural to assume that disorder, in addition to pro- 
ducing noise, must introduce a quasiparticle scattering 
rate 7 which broadens <5p(q, lj), but it is in fact not ob- 
vious. Formally, the elastic scattering rate appears when 
one averages a large number of impurity configurations, 
and is not defined when one studies a particular config- 
uration of impurities. If, for example, we try to take the 
disorder average (p(q, lu)), where (...) indicates a con- 
figuration average, we find that since /(q) is a sum of 
random phases, (/(q)) = and 



(p(q,w)) = [Npo(u) + JVify(q,a;)]<$ q ,o. 



(3) 



In other words, one does not simply find that p(q, uj) 
is broadened by disorder, but rather that in the process 
of configuration averaging p(q, uj), one loses the finite-q 
signal altogether. 

It is therefore the goal of this paper to develop a sys- 
tematic disorder-averaged formalism for understanding 
many-impurity effects in FTSTS. In Sec. Ill Al we demon- 
strate that the variance, (Ap 2 (q, w)) = (|p(q, w)| 2 ) — 
|(p(q, oj)}\ 2 , retains the useful finite-q information when 
disorder-averaged, and that the random noise described 
by /(q) is eliminated as a result of configuration averag- 
ing. Results are illustrated with numerical calculations 
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FIG. 1: (color online) Effect of randomness in the impurity 
distribution on the FTDOS of the 2D tight-binding metal de- 
scribed in Sec. Ill Bl (a) p(q, u>) is plotted for a single pointlike 
impurity embedded in a 100 x 100 lattice. The Fermi surface 
contour = is shown in (b). Arrows qi and q2 show 
the scattering transitions producing the features labelled qi 
and q 2 in (a). The FTDOS for 3% randomly distributed 
impurities is shown averaged over iV c f g configurations for (c) 
jVcfg = 1, (d) 7V cfg = 4, (e) 7V cfg = 10, (f) JV cfg = 29. In (c), 
|p(q,tj)j is plotted, while in (d)-(f), we plot the root-mean- 
square FTDOS p rms (q,o;) = J (Ap 2 (q,oj)} jv cfg . In all figures, 

we take ui — and an impurity potential of V = 10ti. Other 
model details are given in 

Sec. EH 



described in Sec. Ill Bl Finally, we show in Sec. Ill Cl that 
one can develop a diagrammatic perturbation theory for 
(Ap 2 (q, u>)) in which many- impurity effects can be un- 
derstood systematically. One useful consequence of these 
calculations is a demonstration that the mean free path 
can be extracted from FTSTS spectra. A summary of 
these results will be presented in Sec. IIIII 
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II. THEORY AND CALCULATIONS 

A. Simple Model 

For a finite number -/V c f g of impurity configura- 
tions, the variance of the FTDOS is (Ap 2 (q, oj))n f = 
fg _ w ))iVcfg| 2 - We keep N c{g finite for 
the moment because, in practical applications where it 
can be time consuming to obtain spectra of the disor- 
dered system, one might find that it is only possible to 
average over a relatively small number of configurations, 
and it is therefore useful to understand how the noise 
scales with N c f s . Configuration averaging is performed 
by averaging over sets of impurity positions denoted by 
R™ with n = 1, . . . , iVcfg indicating the different configu- 
rations, and i = 1, . . . , Ni labelling impurities. We typi- 
cally need to compute averages of the form 



E< 



-tqR, 



, Wcfg Ni 
& n— 1 2—1 




FIG. 2: (color online) Scaling of FTDOS with number of im- 
purities. Cuts of prms(q, ui) / \TN~i along q = (g, 0) are shown 
for different impurity concentrations and for iY c f g = 29, uo — 0. 



Since, for fixed q, this is just a random walk in the com- 
plex plane with N c f g Ni steps, it follows that 



/ Ni 

£' 



N* s 



N, 



cfg 



/(q), (4) 



where, as always, /(q) is a random complex function of 
q for q ^ with |/(q)| ~ 0(1) and /(0) = 0. Using 
Eqs. (0 and J1J, one can show that 

(p(q,w))iv cfg = [Np (w) +NiSp(q,uj)]Sq }0 



+5p(q,u)* 



cfg 



/(q), (5) 



and that 



(|p(q,^)| 2 )iv ofg 



[iVp (w)+iV i (5/j(q,w)] 2 5 q> o 
Re/(q) 



1 



cfg 



(6) 



B. Numerical Method 

In order to illustrate the results of the previous section, 
and to compare with the diagrammatic perturbation se- 
ries presented in the next section, we have calculated 
the FTDOS exactly for noninteracting quasiparticles in 
a disordered 2D metal. For simplicity, we take pointlike 
disorder and a tight-binding band structure of the form 



Ck 



to — 2t\ (cos k x + cos k y ) + 4^2 cos k x cos k y 



This band structure is often used to describe high tem- 
perature superconductors, but is chosen primarily be- 
cause it is amenable to the real-space numerical method 
described here. For concreteness, we measure energies in 
units of t\ and take to = 0.5ii and t% — 0.45ti. The Fermi 
surface for this band structure is shown in Fig. ^b). 

In real space, the Hamiltonian is defined on a square 
lattice 



Ni 



^ = EE' 



V C ria c rjCr 



(8) 



c=i 



so that 



(Ap 2 (q,c)> 



N a - 



Ni\Sp(q,uj) 



Re/(q) 



cfg 



(7) 



to leading order in -/V c f g . Equation (JJJ is the principal re- 
sult for this subsection, and it demonstrates the key ideas 
of this paper which are that (i) the signal-to-noise ratio 
of the variance of p(q, uS) improves with configuration 
averaging and (ii) limjv cfs ^oo(Ap 2 (q, u))n c{ii is nonvan- 
ishing. The latter point suggests that one can construct 
a disorder-averaged perturbation series for (Ap 2 (q, us)), 
which we will do in Sec. Ill CI 



where c ri<T is the annihilation operator for a quasiparticle 
on site i with spin cr, and ta = to, Uj = —t\ for i and 
j nearest neighbor sites, and Uj = t<i for next-nearest- 
neighbor sites. As before, R^ are the impurity site loca- 
tions, and the impurities have strength Vo = lOij.. Equa- 
tion (JSJ) can be diagonalized numerically, and the eigenen- 
ergies E n and eigenvectors ^(r^) can be inserted into 
Eq. (JTJ to find p(r, u), from which Ap 2 (q, ui) follows di- 
rectly. We remark that for simulations in which N c f s = 1 , 
we plot |p(q,w)|, while for iV c f g > 1, we actually plot the 
root-mean-square variance 



p rms (q,w) = sj (A i o 2 (q,w)) A r cl 
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which is most the most natural quantity to compare with 

|p(q.^)l- 

The numerical approach adopted here restricts the sys- 
tem size to < 10000 lattice sites, and other methods^ ex- 
ist which allow one to study larger systems in the dilute 
impurity limit. However, in many interesting situations 
(eg. ones in which interactions are important) numeri- 
cal calculations will be similarly restricted. This system 
size is also comparable to typical surface areas studied 
in FTSTS experiments. In general, one can expect the 
FTDOS to be consistent with the predictions of the per- 
turbation theory introduced below provided the system 
is larger than the mean-free path. 

The scaling with N and -/V cfg suggested by Eq. has 
been tested numerically using the method outlined in this 
section and holds for a wide range of impurity concentra- 
tions, as shown in Fig. ^ and Fig.|2| In Fig.^ one sees 
that an average over only four configurations improves 
the signal-to-noise ratio significantly (by about a factor 
of 2) and that fine details of the FTDOS can be resolved 
when one averages over w 10 configurations. One also 
sees in Fig. |2 that the spectrum generally scales as y/Nl 
as expected. Closer inspection reveals that most of the 
peaks evident at smaller impurity concentrations are due 
to finite size effects and disappear at larger concentra- 
tions where the mean-free-path is less than the system 
size. Furthermore, the peak at q w (1.9,0), correspond- 
ing to qi in Fig. ^ deviates from the simple scaling at 
smaller values of Ni, but satisfies the ^/N-scaling rela- 
tion best when the mean-free-path is less than the system 
size. 



C. Diagrammatic Method 

The goal of this section is to develop a diagrammatic 
perturbation series for (A/? 2 (q, u>)). The LDOS is defined 
in terms of single-particle Green's functions as 

p{r,uj) = -— [g + (r,r,u>) - 0"(r,r,cj)] , 

where Q + and Q~ are the retarded and advanced Green's 
functions for a particular disorder configuration. Then we 
can write, 

l^-)| 2 = 7^2 E *i S2 £ e ~ 4q - ri2 

si,s 2 =± ri,r 2 



v E **** E 



(2tt) 

si,S2=± ki,k 2 
x^ 1 (ki, kx - q, u)g s - (k 2 , k 2 + q, w) 

where ria = ri — r 2 . The disorder average of is 

(g ± (k,k + q,iu)) = G ± (k,u})S n , , 



kj-q 



5 q,0 



(a) 



(b) k, 



(T) (T) 



(c) 
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FIG. 3: Diagrams relevant to the calculation of (Ap 2 (q, w)). 
(a) Disconnected diagrams cancel in the difference 
<|p(q,w)| 2 ) - |(p(q,a;))| 2 . (b) Leading order and (c) 
second order diagrams which are computed in the text. 
In all diagrams, solid lines represent fully self-consistent 
Green's functions (retarded or advanced) and each pair of 
dashed lines connecting to a vertex labelled "T" represents 



where G(k, us) is the usual self-consistent disorder- 
averaged Green's function. It follows that 

2 



<|p(q^)l 2 ) = 



s,k 



Sq,0 



1 



(2tt) 2 



E siS2 E 



si,S2=± ki,k 2 

x((^(ki,k! - q,c^ 2 (k 2 ,k 2 + q 3 w))), 

and ((. . .)) refers to a disorder-averaging process in which 
only connected diagrams are retained. The first term 
on the right hand side of this equation corresponds to 
the disconnected diagram shown in Fig. |3Ja), and is also 
equal to |(/o(q, oj))\ 2 . Thus, the variance of p(q, w) is 

<Ap 2 (q,c)) = --^ £ SlS2 J^(q,u,) (9) 

where the disorder-averaged two-particle response kernel 
J slS2 (q,w) is 

1 = m E « gsi ( k i- k i - * ^ gS2 ^ k2 + * 



ki,k 2 



(10) 

We first consider the leading order contribution to the 
FTDOS, which is shown in Fig. EJb) and which can be 
written as 



(Ap 2 (q,^)) = Ni 



Im 



^E G+ ( k »-)^ k+q H 



xG + (k + q,w) 



(11) 



where T^ k+Cl (ui) is the retarded scattering T-matrix de- 
scribing scattering from a single impurity embedded in a 
medium with impurity concentration rii = Ni/N. The 
astute reader will recognize that this appears to sim- 
ply be the square of single-impurity result (multiplied 
by Ni) which is widely used to describe the FTDOS 
in strongly-correlated superconductors^ The key differ- 
ence here is that the Green's functions and T-matrix 
are now evaluated at the level of the self-consistent T- 
matrix approximation 16 (SCTMA). This approximation 
amounts to solving the one-impurity problem in an effec- 
tive medium characterized by a self-energy E(k, uS) which 
approximately describes the remaining Ni — 1 impuri- 
ties. In this approximation, the self-consistency condi- 
tion is S(k, uj) — niTk t k(u>), and it can be shown that 
this is equivalent to the coherent potential approxima- 
tion in the dilute limit. This result — that the addition of 
a self-energy to the Green's function accounts for many- 
impurity effects to leading order — has been taken as an 
ansatz by several authors but has not previously 
been formally justified. 

The evaluation of Eq. (|11[) is simplest in the case of 
a random distribution of pointlike impurities with scat- 
tering potential potential Vq. In this case, the T-matrix 
loses its dependence on k and q and we have the follow- 
ing system of equations which must be self-consistently 
solved: 



E ± (w) = n i T ± (o;) 



G ± (k,w) = [w* -e k -S ± (o;) 

k 

T±H = V /[l -V o9 ± (uj)} 



with uj^ — ui ± ir\ and i] a positive, real, infinitessimal. 
The results of this calculation are compared with numer- 
ical results of the previous section in Fig. In general, 
the agreement between the leading-order diagrammatic 
and numerical approaches is good, except for the broad 
peak centered at q — which is present in the numerical 
calculations and absent in the leading-order approxima- 
tion. This peak is distinct from the q = peak in Eq. i ffl l 
(which occurs strictly at q = and is ~ O(N)) and arises 
from higher order diagrams. The finite-q peak heights in 
the diagrammatic calculations scale as y/Nl, as shown 
in Fig. 0Jc) and (d), in agreement with the scaling of 
the numerical results discussed previously. The scatter- 
ing rate in (d) is approximately twice that in (c), and it 
is interesting that, apart from the overall scaling of the 
spectrum, the most significant effect of scattering is to 
enhance the background. Thus, one's ability to resolve 
the peaks is not limited by a reduction of their height, 
but by a reduction in contrast with the background. 

In order to get an understanding of the effects of higher 
order diagrams, we consider the three possible second 




FIG. 4: (color online) Comparison of diagrammatic and nu- 
merical /9rms(q, u)/v~Ni at u> = 0. Results are shown for 3% 
impurities for (a) the numerical method with iV c f g = 29, and 
(b) the leading order diagrammatic result given by Eq. itTTl . 
Cuts along the line q = (q, 0) are shown in (c) and (d) for 
3% and 6% disorder respectively for exact (red triangles) and 
diagrammatic (blue circles) calculations. 



order diagrams shown in Fig.^c), which combine to give 
N,T S1 T S2 



/| 1S2 (q) 



N 2 



£ G Sl (k 1 )G Sl (k 1 -q) 

ki,k 2 



xG S2 (k 2 )G S2 (k 2 + q)[IF SlS2 (k 1 + k 2 ) 
+W S1S2 (ki - k a - q) + 2W S1S1 (q)] (12) 

W SlS2 (p) = mT Sl T S2 — ^2G Sl (k)G S2 (k + p), (13) 



where we note that Ty ++ (p) = W (p)* and W + (p) = 
W~ + (p). Both W ++ (p) and W + ~(p) are shown in 
Fig. at w = 0. The term W+-{p) « 1 - Dp 2 / 2 7 
for small p, where D is the diffusion coefficient of the 
metalfi^ and is substantially larger than l4 7++ (p). 

We are interested in the largest contribution from the 
second order diagrams, which comes from the retarded- 
advanced (W^ ) channel. Because W A (p) is peaked 
at p = 0, we can evaluate Eq. © approximately as: 



7T N 



E|G + (k)G+(k + q)| : 



(14) 



The results of this calculation are shown in Fig. Etc). 
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FIG. 5: (color online) Second order corrections to p rms (q, uj). 
Panels (a) and (b) show the vertex functions W~* (q, lo) and 
W ++ {<\, uj) respectively from which second order corrections 
are calculated, (c) shows p rms (q,tj) evaluated using Eq. 114ft : 
(d) shows a least-squares fit of Eq. 116|l to the results of exact 
numerical calculations at uj = 0; (e) shows the scattering rate 
7 extracted from the SCTMA self-energy (square) and from 
fits of Eq. H6H (solid circle) to the q = peak of numerical 
calculations. Error bars are the standard deviation between 
the fitted curve and data. Panels (a)-(c) are at lo — and for 
m — 0.03, while (d) and (e) are for m — 0.06 and JV c f g = 29. 



One sees a weak background structure which is essentially 
the same as that from the leading order contribution in 
Eq. but the dominant feature in the second order 

correction is the central peak. Physically, this suggests 
that two-impurity processes generally interfere destruc- 
tively and tend to shift weight to small values of q. At 
this level of approximation, however, we do not capture 
the physics of trapped states, which have been predicted 
to occur in d-wave superconductors between pairs of im- 
purities satisfying an appropriate resonance criterion. 20 
Such resonances would tend to enhance finite-q peaks in 
the FTDOS. 

We note that, while Fig. [SJc) has the same qualita- 
tive structure as the exact numerical results shown in 
Fig. Q, the quantitative agreement is still poor. This is 
because, in the absence of inelastic dephasing or coupling 
to the third dimension, weak localization physics affects 
the FTDOS. One can attempt to go beyond second order 
in the perturbation series, and the natural step to take is 
to sum either the ladder or maximally crossed diagrams. 



Both of these calculations give the same result: 
(Ap 2 ( q ,^))~£-^]r|G+(k)G+(k + q )| 2 . 

P P k 

In 2D, the sum over p diverges logarithmically, but the 
divergence will be cut off either by inelastic dephasing or 
by tunneling out of the surface states^ What is perhaps 
most interesting about this result, however, is the q-space 
structure of the central peak, which is the same as in 
Eq. |H|). 

We can estimate the form of the central peak in 
(A,o 2 (q, u>)) in the limit q 2 /2m <C vpq, where for an 
isotropic Fermi surface: 



l^|G + (k)G+(k + q)| 2 = 

k 

such that, for small q, 

(Ap 2 (q,uj)) = 



irN 



7 2 V(^9) 2 + (27) 2 ' 
(15) 



A 



(16) 



where I — vfM is the mean- free path and A dictates 
the peak weight. This result suggests that one can ex- 
tract I from FTSTS experiments. As a check, we have 
made least-squares fits of Eq. Ijlfijl to the central peak 
of an exact numerical calculation of the type outlined 
in Sec. iHBl with m = 0.06 and N cig = 29. We have 
treated both A and £ as fitting parameters, and an ex- 
ample of the numerical data and fitted curve shown in 
Fig. Eld) illustrates that Eq. H16|) appears to describe 
the exact result rather well. In Fig-Efe), the scattering 
rate 7 = vf/( extracted from the fitting procedure is 
shown as a function of uj. For comparison the scatter- 
ing rate 7 = — Im Y, + (uj) calculated using the SCTMA 
is also shown. Because the band structure is anisotropic, 
we have chosen to fit p(q, uj) as a function of the radial 
coordinate q — |q| (note that the point at q = must 
be excluded), as shown in Fig. |SJd). The fit is therefore 
to an angular average of the data, and the error bars in 
Fig. Ef e) are the root-mean-square deviation between the 
fit and the data. We have experimented with fitting to 
cuts through the data along various directions, for exam- 
ple q = (q, 0) and q = (q 7 q), but have generally found 
this to be less reliable than fitting to the angular aver- 
age. In general, since we are attempting to fit a peak of 
half- width w 2£~ 1 , a reliable fit requires a q-space reso- 
lution Aq < 2£ _1 , or (equivalently) that the region over 
which one makes the LDOS map is larger than roughly 
twice the mean-free-path. With these caveats in mind, 
the good agreement between the self-consistent T-matrix 
and exact numerical results is heartening, and shows that 
one can legitimately hope to extract the mean-free path 
from experiments. 

Throughout this work, we have presumed that one is 
averaging distinct samples. However, for experimental 
purposes, it may be more straightforward to measure a 
single large surface area (of side L), and partition the 



7 



data into n 2 subregions of side L/n. The benefits of this 
procedure are that one gains a factor of n in the signal- 
to-noise ratio, but the cost is a loss of q-space resolution 
by the same factor n. Provided L/n > £, however, one 
is not losing important information as a result of the 
partitioning process. 

We finish this section with a brief discussion of inelastic 
scattering. The essential point is that, since the correla- 
tion function measured by Ap 2 (q, uu) is between single- 
particle Green's functions, and is not a two-particle cor- 
relation function, there are no "vertex-correction" dia- 
grams in which inelastic propagators connect the upper 
and lower Green's functions in Fig. [3] To put it an- 
other way, unless one has localized bosonic modes at- 
tached to the impurity sites, inelastic scattering will be 
spatially uniform and will not produce Friedel oscilla- 
tions. Such Friedel oscillations would be described by 
inelastic vertex-correction diagrams, and thus these di- 
agrams must be absent. The effect of inelastic scatter- 
ing is, therefore, solely to modify the single-particle self- 
energy. Since it is the total scattering rate which deter- 
mines peak heights and widths in Ap 2 (q, oj), separation 



of elastic and inelastic scattering rates is difficult. In 
particular, the scattering rate one extracts from Eq. i|ltj|) 
contains both the elastic and inelastic scattering rates. 



III. SUMMARY 

We have studied the disorder-averaged variance 
(Ap 2 (q, u>)) of the density of states of a disordered 2D 
metallic surface. Through a combination of numerical 
and analytic calculations, we have shown that disorder- 
averaging is an effective means of increasing the signal- 
to-noise ratio in the Fourier transformed local density of 
states. This resolves a problem with the conventional ex- 
perimental technique of measuring p(q, to) for which the 
signal-to-noise ratio remains fixed as one averages over 
impurity configurations. Furthermore, we have shown 
that one can construct a diagrammatic perturbation the- 
ory for (Ap 2 (q, oj)) and that it is possible to extract an 
important physical property, the mean-free-path, from 
the central peak in the FTSTS spectrum. 
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